Front motion in an A + B ^ C type reaction-diffusion process: 
Effects of an electric field 
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We study the effects of an external electric field on both the motion of the reaction zone and the 
spatial distribution of the reaction product, C, in an irreversible A~ + — » C reaction-diffusion 
process. The electrolytes A = {A'^,A~) and B = {B^,B~) are initially separated in space and the 
ion-dynamics is described by reaction-diffusion equations obeying local electroneutrality. Without 
an electric field, the reaction zone moves diffusively leaving behind a constant concentration of C-s. 
In the presence of an electric field which drives the reagents towards the reaction zone, we find that 
the reaction zone still moves diffusively but with a diffusion coefficient which slightly decreases with 
increasing field. The important electric field effect is that the concentration of C-s is no longer 
constant but increases linearly in the direction of the motion of the front. The case of an electric 
field of reversed polarity is also discussed and it is found that the motion of the front has a diffusive, 
as well as a drift component. The concentration of C-s decreases in the direction of the motion of 
the front, up to the complete extinction of the reaction. Possible applications of the above results 
to the understanding of the formation of Liesegang patterns in an electric field is briefiy outlined. 

PACS numbers: 05.60.Gg, 64.60.Ht, 75.10.Jm, 72.25.-b 
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I. INTRODUCTION 

Nonequilibrium systems exhibiting pattern formation 
are ubiquitous in nature and these patterns often emerge 
in the wake of a moving reaction front jlj. A classical 
example which motivated our present study is the emer- 
gence of Liesegang bands [^|^ associated with the generic 
reaction scheme A + B —> C and the subsequent precipi- 
tation of C-s. In a typical experimental setup, a chemical 
reactant, B (called inner electrolyte), is disolved in a gel 
matrix, while a second reactant, A (outer electrolyte), of 
much higher concentration is brought into contact with 
the gel. The outer electrolyte diffuses into the gel, reacts 
with the inner electrolyte, and the reaction front moves 
into the gel. Under appropriate conditions, the reaction 
product precipitates and, depending on the geometry of 
the experimental setup, one can observe the emergence of 
families of bands or rings perpendicular to the direction 
of motion of the front 3, 4J. 

This pattern formation process is a rather complex one, 
due to the delicate interplay between the motion of the 
reaction front and the precipitation dynamics of the re- 
action product C. Its study has a history of more than 
a century and the generic, empirical laws describing 
it are well established j^. There is still disagreement, 
however, about the mechanisms underlying this pattern- 
forming process (see |^ Q for a comparative analysis of 
the existing theories). Moreover, with the exception of 
the recently-proposed phase separation (or spinodal de- 
composition) scenario .S, :.9j , the existing theories are es- 
sentially qualitative, as they contain parameters that can- 
not be inferred from experiments, and therefore they do 
not allow for quantitative predictions or for any direct 
comparison with the existing experimental data. Having 
this state of affairs, we believe that it is important to 



test the spinodal decomposition theory in as many ways 
as possible. 

The basic ingredients of the spinodal decomposition 
approach are the presence of a moving reaction front 
and the phase separation that takes place in the wake 
of the front. The dynamics of the inert reaction prod- 
uct C is assumed to have no feedback on the dynamics 
of the reagents ^^l; thus one represents the dynamics 
as a two-stage process. Namely, (i) the production of 
C in the moving reaction front is described by reaction- 
diffusion equations, while (ii) the diffusion and precipita- 
tion of C particles is modelled as phase separation (spin- 
odal jiecomposition) llj using the Cahn-Hilliard equa- 
tion [l^ |lj| , with a source term corresponding to the 
production of C by stage (i). All the parameters involved 
in this model (including the parameters of the Landau- 
Ginzburg free energy associated with the dynamics of C 
particles) are either directly accessible experimentally or 
can be inferred from experimental data (see (3). 

The first stage {A + B ^ C) oi the process appears 
in many other physical and chemical processes and has 
been much studied. The dynamics of the front and the 
spatial distribution of the rate of the production of C (the 
source for the second stage) are known 15] for the case of 
neutral reagents A and B. It should be realized, however, 
that A and B are usually electrolytes which dissociate. 



A^A+ + A- 



B ^ B^ 



B- 



(1) 



and the basic reaction process takes place between the 
'active' ions, e.g.. 



A- 



C, 



(2) 



while the 'background' ions A'^ and B~ are not reacting. 
For the first sight, the background ions may be not im- 
portant. They just ensure local electroneutrality and the 
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Debye screening length is much smaher than the spatial 
lengths involved in the macroscopic pattern formation. 
However, as it was shown recently the dynamics of 
the background ions may generate macroscopic effects 
(even if the screening length is negligible) and may influ- 
ence both the propagation of the reaction front and the 
structure of the resulting pattern. Thus, an obvious way 
to obtain more insight into the details of the dynamics 
of Liesegang pattern formation (with the ultimate goal 
of validating one theoretical scenario or another) is to 
switch on an electric field and to study its effects. 

Experimental and theoretical investigations in this di- 
rection have been going on f or q uite a while, see [TtLITI . 

El ES IM 111 S Im il 111 Ea El 111 113 . The emerg- 
ing picture, however, might be sometimes confusing. For 
example, on the experimental side. Das et al.'s experi- 
ments in 113 show a diffusive motion of the front for 
a polarity of the applied field that favors the reaction, 
and for a reverse polarity, as shown by Sultan et al. |27j 
the motion of the front acquires a supplementary drift 
component. These results are just the opposite of those 
obtained by Lagzi 29\, however it can be argued 31] that 
in the latter case the properties of the intermediate com- 
pounds are responsible for this 'anomalous' behavior. As 
far as the theories are concerned, a drift of the front is 
obtained analytically in e.g. |13 under the irrealistic as- 
sumption of the constancy of the electric field along the 
system, or it is put 'by hand' as e.g. in [soj . 

In order to clarify the effect of the electric field on 
Liesegang patterns, we propose to revisit the problem in 
the framework of the spinodal decomposition theory. In 
the present paper, we shall concentrate on the first stage 
of this scenario, namely on the effect of the electric field 
on the motion of the reaction front and on the concentra- 
tion of the reaction product C left behind the front. As 
mentioned above, the problem of the first stage is more 
general and thus we feel that the results are of interest in 
their own. Although we shall mention the implications 
of the reaction zone results for Liesegang patterns, the 
calculations and discussions of the second stage of the 
process will be presented in a forthcoming study. 

As it will be discussed in the next section (Sec. m 
with some details in Appendix A), the new element 
in our model is that the dynamics of all the ions 
(yl+, , B^, B^) is followed by using reaction-diffusion 
equations obeying the electroneutrality condition. Fur- 
thermore, the electric field is taken into account by pre- 
scribing a given potential drop between the two ends of 
the sample, as in real experimental situations. The equa- 
tions are solved numerically in Sec. IIIII where, in view of 
the possible comparison of our results with experimen- 
tal data, we concentrate our study on the experimentally 
relevant range of parameters and observation times. 

Our results can be summarized by contrasting them 
with the fieldless case where the reaction front moves 
diffusively, and it leaves behind a constant concentration 
Co of the reaction product C 0, 0|. In the presence 
of the electric field, the polarity of the field is an im- 



portant factor. If the field favours the reaction (in the 
sense that it drives the reacting ions A~ and B^ towards 
the reaction zone) , then the motion of the reaction front 
remains diffusive, and the concentration of the reaction 
product c{x) increases linearly in space in the direction 
of the motion of the front. The diffusive front motion is 
in agreement with Das et al.'s experiments HI E^ and 
in contrast with theoretical conjectures about constant- 
velocity drifts m Ell . An important remark is that this 
case is difficult to handle, since a strong electric field de- 
velops in the reaction zone (especially for large applied 
tensions), which leads to numerical instabilities and thus 
restricts our numerical solutions to relatively short times. 

In the case of opposite polarity of the applied field, 
the motion of the reaction zone is still diffusive at short 
times, but a crossover to drift at large times, as seen in 
Sultan's experiments |27j , can be inferred from numerical 
data. As far as the concentration of the reaction product 
is concerned, we find that at small times it decreases 
linearly in the direction of the motion of the front, while 
a slow-down of this decrease is noticed for larger times, 
up to the extinction of the reaction. 

The conclusions and the implications of the results for 
Liesegang phenomena will be presented in Sec. IV. Fi- 
nally, some salient details of the problem such as the 
electric-field profiles, the origin of the numerical insta- 
bilities, the time evolution of the electric current, the 
finite-size effects, etc. are discussed in Appendix B. 



II. THE MODEL IN THE PRESENCE OF AN 
APPLIED ELECTRIC FIELD 

The experimental setup we shall consider is repre- 
sented schematically in Fig. ^ The two electrolytes are 
initially separated from each other, each is uniformly 
spread into a gel matrix, and the concentration of A's 
is much larger than that of -B's. Thus the reaction front 
A~ + B^ C moves to the right along the column. An 
appropriate choice of the experimental conditions (type 
of reagents, concentrations, etc.) leads to quasiperiodic 
precipitation of the inert reaction product C in the wake 
of the moving front - the Liesegang patterns, see Fig. ^ 
Due to the presence of the gel the convection phenomena 
are absent, and thus the evolution of the concentration 
fields of the particles is accurately described by reaction- 
diffusion equations. 

Our model is based on several simplifying assumptions, 
and we refer the reader to Appendix A for more details 
on a more general model. 

(i) The first assumption is that the system is one- 
dimensional, which means that all the relevant quantities 
depend on a single spatial coordinate x, with —La < x < 
Lb being the spatial extent of the system. This is essen- 
tially realized in practice provided that the length of the 
gel column is much larger than its width. The edge ef- 
fects in the transverse directions rnight have some slight 
consequences (e.g., as reported in |22| . in the presence of 
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FIG. 1: Schematic representation of the system under study. 
The electrolytes A (^+, A~) and B (_B+, B~) are initially lo- 
cated in the regions [~La, 0), respectively (0, Lb]- The pre- 
cipitation bands, i.e., the alternation of high-density-C re- 
gions (shaded areas) and low-density-C ones, emerge in the 
wake of the moving reaction front (dashed line at Xf{t)). For 
an applied tension U, an electric current of density is 
flowing through the system. 



an external electric field the diameter of the tube has a 
small influence on the propagation of the reaction front). 
However, the final pattern is still one-dimensional to a 
very good accuracy. 

(ii) We consider instantaneous 100% dissociation of the 
electrolytes A and B according to Eq. |^ (the assump- 
tion of "ideally strong" acid and basis). Correspondingly, 
only and ions (and the reaction product C) are 
present in the system. 

(iii) The ions A~ and B^ are reacting irreversibly with 
an infinite reaction rate. This can be justified by the 
fact that the characteristic time of the reaction process 
is 3 to 6 orders of magnitude less than the time scales 
invloved in the diffusion and precipitation processes, see, 
6-g-, 0, ll^ . In this case, one has a pointlike reaction 
front at a time-dependent position x j [t) |3a|. 

(iv) The electroneutrality approximation p33 | - accord- 
ing to which the local charge density is zero on space 
scales that are relevant to pattern formation - reads 



Zi ni{x, i) = , 



(3) 



where zi is the charge of the i-ih ion - in terms of the 
elementary charge q - and nAx,t) denotes its density. 
As discussed in detail in [T^ 13^ , this condition is well 
fulfilled in the systems we are investigating. Moreover, 
we shall consider monovalent ions, \zi\ = 1 for all i. 

(v) Finally, we consider equal diffusion coefficients of 
the ions, Di = D for all i. 

The initial conditions for the concentration profiles 
a^{x,t) and b^{x,t) of the A^ and B^ ions correspond 



a {x,t = 0) = a+(x, t = Q) — ao , 
b+{x,t^ 0) = b-{x,t = 0) = 0, 

for — La < X < 
a^{x,t = 0) = a+(.T,t = 0) = 0, 
6+(x,t = 0) = b-{x,t = 0) = 5o, 

for 0<x <Lb, (4) 

with flo 3> bo. We also suppose that the concentrations 
of the ions arc maintained fixed at the borders, 

a+(-L^,t) = a^{-LA,t) = ao , 
b+{-LA,t) = b-i-LA,t)^0, 
a+{LB,t) ^ a-{LB,t) ^ 0, 
b+{LB,t) = b-{LB,t) = bo, (5) 

i.e., the system is in contact at its left and right borders 
with two 'infinite' reservoirs of ions A^ and B^, respec- 
tively. 

Under these assumptions, the general evolution equa- 
tions for the concentration profiles a^{x,t) and b^{x,t) 
of the A^ and B^ ions presented in the Appendix A 
acquire the following simplified expressions: 

(i) For X < Xf{t) (i.e., in the wake of the reaction front) 



da {x,t) 
di 

6+(.T,<) = 
da+{x,t) _ 
dt ~ 
db-ix,t) 



D 



dx^ 



D 



= D- 



dx^ 



2q dx \ a+ 



ill. 

2q dx \ a^ 



dt dx^ 
(ii) For X > Xf{t) (i.e., ahead the reaction front) 



a {x,t) = , 
db+{x,t) 
di 

da~^{x,t) 

di 

db-ix,t) 



D 



d^b+ Jo d fb+ 



D 



D 



dx^ 
d^a+ 
dx^ 

d%- 



2q dx \b' 
jo_d_ f a_ 
' 2qdx\ b- 



(6) 
(7) 
(8) 

(9) 

(10) 
(11) 
(12) 

(13) 



dt dx^ 
Here jo is the time-dependent electric current density 

-2qFDU 



Jo{t) 



RT 



dx 



(14) 



xfW 



determined by the constant voltage difference U = 
V{Lb) — V{—La) apphed between the two ends of the 
system 35]. F = qNA is Faraday's constant (i.e., the 
electric charge transported by a mole of monovalent pos- 
itive ions), R is the universal gas constant, while T is the 
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temperature. Finally, the local electric field is given by: 
Jo{t) 



E{x,t) 



qFD 
RT 



i 

-2U 



E 



r^'*^ dx 


/■^« dx\ 




Jxfit) b- J 



(15) 



At a first sight, one would be tempted to conclude 
that the equations (|S|)- (|13|I for the concentration fields 
look simple, but actually (i) they are nonlinearly cou- 
pled through the integral quantity joit), according to 
Eq. 1141) . and (ii) they acquire time- dependent boundary 
conditions, 



a-{xf{t))=b+ixf{t)) = 0, 



(16) 



which are complicated to handle. The meaning of these 
latter conditions is just that the concentrations of the 
reagents are zero at the front, and their flux towards the 
reaction front are equal. 

Even under the simplifying assumption presented 
above, an analytic solution for the concentration profiles 
of the ions and thus for the motion of the front seems 
hopeless. We thus turned to numerical solutions, whose 
main results are presented in the following Section. 



III. NUMERICAL RESULTS 

As already mentioned in the Introduction, we decided 
to focus our analysis on situations that are experimen- 
tally relevant. In particular, although we are aware of 
different other regimes that might appear, we restrict the 
presentation to the ones that are most likely to appear in 
experiments. The choice of the parameters intends thus 
to mimic real experimental situations, namely, we con- 
sidered concentrations of the ions oq and bo of the order 
10"^ M - 10 M, lengths of the system L = La + Lb 
of some tenths centimeters, and tensions U applied be- 
tween system's edges such that the corresponding field 
intensity U/L varies between ±15 V/m. The common 
diffusion coefficient of the ions was chosen as D — 10"^ 
m^/s for all the calculations. Finally, times of observa- 
tion of the system that vary between a few hours and 
some tenth days. Note that this observation time cor- 
responds to the real experimental time and not to the 
time required by the numerical calculations. As far as 
the actual numerical calculations are concerned, we in- 
tegrated the partial differential equations l|^- H13|) using 
a fourth-order Runge-Kutta integration scheme in time 
and a first-order differenciation procedure in space. 

We shall generally present the results that were ob- 
tained for a polarity [/ > of the applied field that fa- 
vors the reaction, i.e., that drives the A~ and ions 



towards the reaction zone. This is the case that is apri- 
ori more interesting for the formation of the Liesegang 
pattern. However, this situation is numerically difficult, 
since a strong electric field develops rapidly (especially for 
large applied tensions) in the reaction zone, thus leading 
to the appearance of numerical instabilities which limit 
the observation time. For more details on this point see 
Appendix B. 

We also indicate the main results for the opposite po- 
larity of the electric field, for which the production of C 
is slowed down and finally stopped. This is a situation 
that becomes rather rapidly uninteresting from the point 
of view of the pattern formation. 



A. Front propagation 

In the limit of an infinite reaction rate, the reaction 
zone reduces to a point whose position Xf(t) is moving 
to the right in the experimental setup of Fig. 

The main result of the simulations we carried on for the 
case of the field favoring the reaction is that the motion of 
the front is diffusive with an effective diffusion coefficient 



Xf{t)^^/2Dji. 



(17) 



The above statement is valid as far as the border effects 
are negligible (see below) and for the parameter domains 
that are experimentally relevant. This result is in agree- 
ment with the rather accurate experimental results of 
Das et al. 22, 23J, but contradicts some previous the- 
oretical considerations [23, HO] which, however, in our 
opinion, have no justification. 

Figure 121 shows the trajectory of the reaction front for 
various tensions U / L applied to the system. A first re- 
mark would be that the effect of the electric field on 
front's motion is very small. One is also confronted with 
a counterintuitive result (that has also been seen experi- 
mentally 22j _23J ) , namely that the motion of the front is 
very slightly slowed down with increasing field (although 
the polarity of this field was chosen to favor the reac- 
tion). We had to limit the observation times for these 
large tensions to 5 hours, because of the numerical in- 
stabilities already invoked above. Of course, for smaller 
tensions one can go to longer observation times, like e.g., 
in Figs. H and 

A look at the concentration profiles of the ions, see 
Fig. 121 allows to sketch an intuitive explanation of this 
unexpected behavior. Namely, if no reaction takes place, 
in the presence of an electric field with a polarity U > 
the negative ions will be pushed ahead the positive 
ones (and the corresponding symmetric situation for the 
ions), as it can be seen from the first panel of the 
figure. With an infinite-rate reaction present, see the sec- 
ond panel of the figure, the A~ ions are consumed at the 
front and the A~^ ions are forced to move ahead, against 
the electric field (a similar, but symmetric situation oc- 
curs for the B^ ions). Therefore, the motion should be 
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FIG. 2: Position of the reaction front as a function of the 
square root of time for different values of the intensity U jL 
of the applied field. The values of the other parameters are: 
La = 1 cm, Lb = 20 cm, ao = 10 M, and &o = 0.1 M. The 
observation time is limited to 5 hours (see the main text). 



diffusive and has to slow down slightly with increasing 
field. 



duction of C) one can neglect the dependence of the dif- 
fusion coefficient Df on the applied tension and consider 
its fieldless expression given by |1] : 




D/\ _ (ao/bo) - 1 
2D I (ao/bo) + 1 ' 



(18) 



that corresponds to an increase in Df with increasing 
concentrations ratio (ao/6o), see Fig. 0] A qualitative 
comparison of Figs. Inland 01 shows that Df is much more 
sensitive to changes in the concentration ratio (ao/^o) 
than to changes in the apphed field U/L. 
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FIG. 3: Two snapshots of the concentration profiles of the 
ions in the presence of an applied electric field U/L — 0.5 
V/m: (a) without reaction (b) with infinite reaction rate. 
The observation time is of 50 hours and the values of the 
other parameters are La = Lb = 100 cm, ao = 10 M, and 
bo = 0.1 M. 

However, for all the practical purposes (e.g., the pro- 



FIG. 4: Position of the reaction front as a function of the 
square root of time for three values of the concentrations ratio 
(ao/bo). The values of the other parameters are: La = 1 cm. 
Lb = 100 cm, and U/L = 0.5 V/m. The observation time is 
10 days. 

The finite spatial extent of the system does not affect 
the diffusive behavior of the front until the reaction front 
reaches the very vicinity of the borders. Then the front 
slows down and stops, and there is a rapid accumulation 
of C product at the border. However, for the usual ob- 
servation times we are dealing with, this regime of the 
front is not relevant, see Appendix B for more details. 

B. The concentration of the reaction product C 

We come now to the most important results from 
the point of view of the subsequent pattern formation 
through precipitation of C. Namely, the study of the in- 
fluence of the applied electric field on the production of 
C in the wake of the moving front. 

As well-known (see, e.g., [EE; ^3); in the absence of 
an applied electric field the concentration of C left behind 
the reaction front is constant, and its value cq is deter- 
mined by the initial concentrations of the ions aq and bo, 
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and by their diffusion coefRcients. In the particular case 
of equal diffusion coefficients D of the ions, its value is 
given by: 



co-aoKj2D/Df, 



(19) 



where K = {1 + 5o/ao)(20F)-i exp{-Df/2D), and the 
diffusion coefhcient -D/ of the front is given by Eq. H18() . 

When an electric field is applied to the system, the 
concentration of C left behind the moving reaction front 
is no longer constant, but becomes space-dependent, c = 
c{x). 

There are two factors that determine the instantaneous 
quantity of C that is produced in the wake of the front, 
namely: (i) The velocity of the front, if — dxf{t)/dt. 
As discussed above, this one is diffusive, with a diffu- 
sion coefficient Df that is essentially unaffected by the 
externally applied electric field, if — \J Df jit. 

(ii) The instantaneous value of the current of the re- 
acting ions A~ and that arrive at the reaction point, 
|ja- = bf)+ that has a diffusive compo- 
nent, as well as a component proportional to the electric 
current density (i) . 

A first regime, that appears without exception when 
the polarity of the applied field favors the reaction, is a 
linear increase of c with x, with a slope that is propor- 
tional with the apphed tension, see Fig.El Moreover, the 
slope is proportional to cq, the concentration of C in the 
absence of the field, given by Eq. (fT^ . 
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FIG. 5: The density of the reaction product C left behind 
the reaction front for different values of the electric field U/L 
applied to the system. The observation time is t = 10 days. 
The values of the other parameters are: La = 1 cm. Lb = 100 
cm, ao = 10 M, and &o = 0.1 M. 



Thus the numerical results can be summarized in 

cix) ^co[l + {a,U/L)x], (20) 

and for the profiles in Fig.|3 one inferres ac ~ 5V^^. 

This behavior of the production of C can be easily un- 
derstood once we notice that for the relevant observation 
times the current density jo is practically constant (see 
below), and thus \ja- ,b+(^f)\ ^ constant ~ jo and a 
diffusive part ^ l/Vt. This fact, in conjunction with the 
diffusive motion of the front, leads to a density of C in 
the wake of the front 

c{xf) ~ Ijq- i,+ {xf)\/if ^ Vi + const. ^ Xf + const. , 

(21) 

i.e., a linear increase with x. 

The regime of constancy of the current density jo is 
discussed in more details in the Appendix B. It is main- 
tained as long as the total electric resistance of the system 
is practically constant. This means as long as the resistiv- 
ity of the depletion (low ion concentrations) region that 
develops around the moving front remains unimportant 
as compared to the resistance of the rest 'unperturbed' 
part of the system. Of course, the longer the system, the 
longer this regime lasts. For most experimental situa- 
tions, this is indeed the only observable regime. 



C. Reverse polarity of the applied tension 

We present below the main results for a polarity of the 
applied tension U < that does not favour the reaction, 
i.e., tends to push the A~ and B"*" ions far apart. A 
first result is that the motion of the front is no longer 
purely diffusive, but, as can be seen from Fig. El a drift 
component emerges. The curves can be well fit using the 
expression 



Xf{t) = ^/2Dji + Vft, 



(22) 



with, however, the drift part playing a significant contri- 
bution only for sufficiently long times and/or high ten- 
sions, in agreement with experimental data p^ . 

A naive attempt to justify qualitatively this behavior 
of the front appeals to the type of argument we have 
already used in Sec. Ill A for the case of the favorable 
polarity of U (in order to explain, there, the diffusive 
motion of the front). Namely, in the case U < 0, the 
A'^ ions always move ahead the (reacting) A~ 's (and the 
symmetric situation for the ions), as seen in Fig. |7| 
but in this case their motion is in the sense of the applied 
field, and might allow for a diffusive, as well as a drift 
component of front's motion. 

The concentration c of the reaction product decreases 
nonlinearly with x, as seen in Fig.|Sl faster when the ap- 
plied field intensity \U\/L is larger, and up to the com- 
plete stop of the reaction (when obviously Eq. H22I) be- 
comes meaningless). Thus this situation becomes rather 
rapidly uninteresting from the point of view of Liesegang 
pattern formation, as it leads to the disappearance of the 
pattern. 
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FIG. 6: The position of the reaction front as a function of 
the square root of time for different values of the apphed field 
U/L < as compared to the fieldless case. One notices the 
drift component of the motion. The observation time is f = 10 
days. The values of the other parameters are: La = 2 cm, 
Lb = 40 cm, ao = 10 M, and feo = 0.1 M. 



FIG. 8: The density of the reaction product C left behind the 
reaction front for different values of the applied field U/L < 
as compared to the fieldless case. One notices the rapid 
extinction of the reaction with increasing field intensity \U\/L. 
The observation time is t = 10 days. The values of the other 
parameters are: La ~ 2 cm, Ls = 40 cm, ao = 10 M, and 
6o = 0.1 M. 
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FIG. 7: Two snapshots of the concentration profiles of the 
ions in the presence of an applied electric field U/L = —0.5 
V/m: (a) without reaction (b) with infinite reaction rate. The 
observation time is of 50 hours and the values of the other 
parameters are La = Lb = 100 cm, ao = 10 M, and bo =0.1 
M. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this paper we examined the influence of an applied 
electric field on the propagation of the reaction front and 
on the production of C for an A~ + 5+ — > C reaction in 
a one-dimensional geometry that is specific to Liesegang- 
pattern experiments. We deliberately concentrated the 
presentation of our results on the experimentally relevant 
parameters and observation times (although, for other 
ranges of the parameters and longer observation times 
other regimes might come into play). Our main con- 
clusions are that for a polarity of the applied field that 
favours the chemical reaction, (i) the motion of the reac- 
tion front is diffusive, with a diffusion coefficient practi- 
cally unaffected by the applied field, and (ii) the concen- 
tration of the reaction product c{x) left behind the front 
increses linearly with x, with a slope that is proportional 
to the applied tension. 

For the reverse polarity (i.e., that unfavours the reac- 
tion), (i) the motion of the front has a diffusive as well as 
a drift component, and (ii) the concentration c{x) of the 
reaction product left by the front decreases nonlinearly 
with X up to the complete stop of the reaction. 

The above results on the motion of the reaction front 
are in agreement with the experiments in p3 . l23ll2^. b ut 
appear to be at variance with the front motion in |29l |. 

The data we presented are based on a set of approxima- 
tions. Any of them can be relaxed and/or improved, but 
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we have evidences that none of them affects essentiaUy 
our conclusions. 

The next step would be to use the above results to de- 
termine the influence of the applied tension on the char- 
acteristics of the Liesegang patterns. This can be easily 
done in the framework of the spinodal decomposition sce- 
nario. The dynamics of the C particles is still described 
by a Cahn-Hilliard equation, but this time the source 
term is modified (as compared to the fieldless case, see 
0) in agreement with our conclusions (i) and (ii) above. 
A rapid examination of these elements already allows us 
to predict for a favourable polarity U > a decrease 
in the spacing between the bands (the more pronounced 
the higher the applied tension), and, for sufficiently long 
times and/or high enough tensions, the disappearance of 
the pattern. For the reverse polarity U < 0, the pattern 
(if any) will present an increase in the spacing between 
bands as compared to the fieldless case, followed by a 
rather rapid disappearance of the pattern. A detailed 
study of this problem will be given in a further publica- 
tion. 



APPENDIX A: THE GENERAL 
ONE-DIMENSIONAL MODEL 

We present below the general one-dimensional evolu- 
tion equations for the concentration profiles of the par- 
ticipating species when we do not the approximations of 
ideally strong acid and basis, infinite reaction rate, and 
equal diffusion coefficients of the ions (see the main text). 

When relaxing the assumption of ideally strong acid 
and basis for the electrolytes A and B, one has to take 
into account their diffusive motion and their finite disso- 
ciation, i.e., the following dynamics of their concentration 
profiles a{x,t) and b{x,t): 

= - K{Kaa a~a+) , (Af) 



db{x,t) 
dt 



D, 



\{Kbb^b-b+). (A2) 



Here a~ {x,t), a^{x, t), b~^{x, t), and b~{x, t) are the con- 
centrations of the appropriate ions resulting from A and 
B, respectively, and Aa,b are the relaxation constants to 
the respective dissociation equilibria characterized by the 
dissociation constants Ka and Kt- 

As already mentioned, A'^ and B^ are not reacting, 
while the ions A~ and B'^ are reacting irreversibly with 
a certain reaction rate k, that we shall consider finite, 



A- 



B^ 



C. 



(A3) 



We assume that the dynamics of the inert reaction prod- 
uct C has no feedback on the dynamics of the reagents 

The macroscopic evolution equations for the concen- 
trations of the ions in the presence of an electric field 



E{x,t), by considering different diffusion coeficients of 
the ions, read: 

= D-Z^ + X^{Kaa-a'a+)-ka-b+ 



dt 



" dx^ 
n-,- F d{a-E) 
" " RT dx 



(A4) 



db+{x,t) _ d%+ 



dt 



^ dx^ 



+ \b{Kbb - b+b-) - ka-b+ 
F d{b+E) 



" ^ RT dx 



da'^(x,t) .d'^a+ , _ 



dt 



" dx^ 



F d{a+E) 
" RT dx 



db-(x,t) d^b- ^ , , , , 

^ = D;:^^ + Xb{Kbb-b+b-) 



dt 



F djb-E) 
'^^'^ RT dx 



(A5) 



(A6) 



(A7) 



Here are the respective diffusion coefficients of the 
ions, F — qNa is Faraday's constant (i.e., the elec- 
tric charge transported by a mole of monovalent positive 
ions), R is the universal gas constant, while T is the tem- 
perature. The Zi-s are signed integers giving the charge 
of the i-th ion, i = a^, 6^, in units of elementary charge 
q, and we hold the local electroneutrality assumption 



Zini{x, t) — Q 



at any point and at any time {ni{x,t) are the concen- 
trations of the ions). It is useful to recall that indeed 
local electroneutrality assumption is well-justified for 
Liesegang type of experiments: at characteristic ion con- 
centrations ^ 10~^ - 10 M usually present in Liesegang 
experiments, the Debye screening length is of the order 
of ~ 10~^° m - 10~^ m; thus, it is indeed negligible as 
compared to the other length scales present in the system 
- lengths of the precipitation zones of the order ~ 10"'^ 
m, and width of the reaction zone ~ 10~^ m. See 0,01 
for a more detailed discussion of this point. 

As known from textbooks [s^, the local electric field 
E{x, t) in the above equations is determined both by the 
externally applied field and by the local electroneutrality 
condition, and it is given by: 



E{x,t) 



duj 
'' dx 



RT ^ 



(A8) 



jo(^) is the electric current density, fiowing through the 
system. Note that in view of the electroneutrality con- 
dition Jo is divergence-free, i.e., for the one-dimensional 
case it is only time-dependent. 
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Suppose that a constant voltage difference U = 
V{Lb) — V{—Lji) is applied between the two ends of 
the system. Then, according to Eq. I)A8|I . 



joit) 



dxE{x, t) 



ox 



dx- 



-La 



-u . 



(A9) 



Thus the instantaneous value of the current density jo {t) 
is determined by the applied tension and by the instan- 
taneous concentration fields ni(x,t): 



FU 



dxY,DrZ~ 



dx 



(AlO) 

At its turn, jo{t) determines through Eq. ljA8|l the instan- 
taneous electric field E{x,t), that enters the evolution 
equations (|X4l) - (|X7|) . 

One concludes that Eqs. ((XlT) . HA2p . (|X4|) . (|X5|) . IjAep . 
(|A7|) for the concentration profiles, IjASp for the electric 
field, and (jAf 0|) for the current density are all coupled 
in a highly nonlinear, intricate way. In our actual study 
(see Sec. II), we used several simplifications to make these 
equations more tractable. 



APPENDIX B: MORE NUMERICAL RESULTS 

Here we shall present additional results from our nu- 
merical simulations. We did not describe them in the 
main text, in order to make the essential results more 
transparent. However, we think they may put supple- 
mentary light on our conclusions, as well as on some con- 
troversial points of previously published results. 

1. Influence of flnite-size effects on front's motion 

This section intends to strengthen the conclusion pre- 
sented in the main text, namely that the finite-size effects 
do not infiuence the motion of the front as long as this 
one does not 'hit' the right border. This reinforce also 
the idea that the factors responsible of front's motion are 
strictly local. 

In Fig. 1^1 we considered three system sizes and observed 
the long-time behavior: (i) for the smallest system the 
front hits the right border during the simulation time; 
(ii) for the intermediate system the front starts to 'feel' 
the border at the end of the simulations; (iii) in the the 
longest system, the front is unaffected by the borders 
during the simulation time. 

Note however the very long time required to see such 
border effects even for the shortest system (about 10"* 



a 




80 120 

1/2 . 1/2 

t (min ) 



FIG. 9: Position of the reaction front as a function of the 
square root of time for different values of the system length. 
The values of the other parameters are: ao = 10 M, bo =0.1 
M, and applied electric field U/L = 0.5 V/m. The observation 
time is of 500 hours. 



minutes for a system with Lb of only 10 cm). In most of 
the real experiments (for which both the system length 
is bigger - of at least a few tenths centimeters - and the 
observation time is usually smaller - a few days), as well 
as in most of our simulations, this slowing-down regime 
of the front is practically never attained. 



2. The current density jo 

For a constant applied tension U the current density 
jo is a function of time only; it is a global quantity, being 
determined by an integral over the entire system involv- 
ing the instantaneous concentration fields. Consequently, 
jo{t) is sensitive to finite-size effects all along its evolu- 
tion. The different regimes of jo(i) result essentially from 
the interplay between the evolving relative electric resis- 
tances of the depletion zone and that of the rest, 'unper- 
turbed part' of the system. 

The depletion zone is defined as the region where the 
concentrations of the ions are significantly smaller than 
their initial values. This region expands progressively, 
more or less rapidly (depending on the applied tension) 
around the moving reaction front. It can develop, of 
course, up to the size that is allowed by the borders of 
the system. It is a region of higher resistivity than the 
'unperturbed' part of the system. There are two ele- 
ments (with distinct temporal evolution) that determine 
the electric resistance of the depletion zone, namely (i) 
its spatial extent and (ii) its degree of depletion, i.e., the 
concentrations of the ions in the region. Note that the 
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resistivity of the unperturbed part of the system is de- 
termined by its spatial extent, which equals L minus the 
length of the depletion zone. 

At the initial stages, the spatial extent of the depletion 
zone is small as compared to the rest of the system, and 
the concentrations of the ions are not extremely low yet. 
Therefore, its contribution to the electrical resistance of 
the column is negligible; during this stage, the resistance 
of the column is essentially constant and determined by 
the 'unperturbed' part of the system. Consequently, one 
obtains an initial regime of constant current density. This 
can be observed in Fig. ^| from which one also realizes 
that this is actually the regime that is accesible experi- 
mentally, and thus the only of interest to us. 



resistivity of the depletion zone becomes comparable to 
that of the rest of the system. Then a transient regime 
starts ~ seen as an increase in for both systems. Fi- 
nally, the increasing resistivity of the depletion zone be- 
comes dominating in the system, causing a monotonous 
decrease in the current density |jo|- However, the in- 
termediate system does not 'have time' to develop this 
regime further: it hits the border and starts the corre- 
sponding increase in |Jo|, according to the mechanism 
sketched at point (i) above. A longer time simulation 
(not shown here) for the longer system shows a continu- 
ation of this regime, with jo{t) ~ 



3. The electric field 
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The electric field is not constant along the system, as 
can be seen from Fig. llH _^in contradiction with the funda- 
mental assumption in 20] . This can be easily understood 
by inspection of Eq. IjASp . taking into account that the 
concentration fields vary along the system. For a polar- 
ity of the externally applied tension U that favors the 
reaction, the absolute value of the electric field increases 
rapidly in the vicinity of the moving front. This effect 
is due to the decrease in both the the concentrations of 
the reacting ions and i?+ (due to the reaction), and 
the concentrations of the background ions and B~ 
(because the polarity of the electric field pushes them 
away from the reaction zone). As a result, the sum in 
the denominator of the expression (jASp determining the 



10000 electric field J2i Dizfui becomes small in the vicinity of 



the reaction front, and thus the local field becomes very 
large. 



FIG. 10: The absolute value of the current density |jo| as a 
function of time (logarithmic scale) for various system sizes - 
the same systems as those in Fig. |^ The values of the other 
parameters are: ao = 10 M, bo = 0.1 M, and applied electric 
field U/L — 0.5 V/m. The observation time is of 500 hours. 



Although we shall not concentrate on them further (as 
not relevant for the typical experimental domains) let us 
briefly describe the other regimes that can be seen in 
Fig. ^1 in the final part of the simulations. 

(i) For the shortest system, the depletion zone can- 
not develop significantly before the front reaches the 
right border. Then the current passes directly from the 
constant-value regime to a rapid increase. When the 
front is in the vicinity of the border, the right-hand (with 
respect to the front) high-resistivity region of small con- 
centration of the ions reduces progressively in size, while 
there is an advancement, behind the front, of a region of 
low-resistivity, high ion concentrations. This generates 
a decrease in the resistance of the system, and thus the 
observed increase in the current density. 

(ii) For the intermediate and the longest system, the 
depletion region has more time to develop. At a certain 
time (shorter for the intermediate system), the electric 
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FIG. 11: Snapshots of the electric field E{x, t) at two times 
and for two different values of the electric field U/L applied 
to the system: (a) U/L = 5 V/m (b) U/L = 15 V/m. The 
values of the other parameters are: La ~ 1 cm. Lb = 20 cm, 
ao = 10 M, and bo = 0.1 M. 
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This rapid - both spatial and temporal, see Fig. 1111 - 
increase in the electric field renders the numerical pro- 
cedure unstable after some time. The other element 
that contributes to this instability is also connected with 
the small densities of the ions around the reaction zone: 
namely, from one integration step to another, the lo- 
cal concentration profiles may acquire unphysical neg- 
ative values. These instabilities develop faster for larger 
applied tensions; their appearance can be delayed by 
decreasing both the spatial and temporal discretization 
steps. One is thus limited in the choice of the upper value 
of the applied tension J7 by a 'reasonable' choice of the 
spatial and temporal discretization steps, that allow the 
study of the system for time intervals that are physically 
relevant. From a technical point of view, these instabil- 
ities are promptly notified by the value of the electric 
field at the reaction point, that acquires an unphysical 
temporal evolution at the moment of the onset of the 
instabilities. 

The appearance of these numerical instabilities can be 
also postponed by: (i) giving up the infinite dissocia- 
tion rate approximation, and allowing the A and B elec- 
trolytes to diffuse towards the reaction zone and to sup- 
ply it with ions and through dissociation; (ii) 
taking into account the finiteness of the reaction rate 
k between A~ and 5+, that leads to a slower decrease 
of the concentration of the reacting ions in the reaction 
zone. Note that in this case the reaction zone acquires 
a finite spatial extent (i.e., it is no longer reduced to a 
single point). 

Finally, we should mention in connection with the elec- 
tric field that, in real systems, this one has supplementary 
inhomogeneities related to the geometry of the system 
and, furthermore, the electric stability of the gel column 
may also reduced by electric effects at the walls of the 
container [3l| . 



4. The concentrations of the ions 

The concentration profiles ni{x, t) of the ions have a 
very complex evolution, that is determined by their diffu- 
sion coefficients, the initial concentrations, the length of 
the system, and the applied electric tension. A system- 
atic presentation of all the corresponding effects would 
be tedious, and thus we shall present below only a few 
relevant elements. 

For a given system (i.e., with fixed parameters), a qual- 
itative typical temporal evolution of the concentration 
profiles corresponds to a reduction of the concentration 
of the ions in the vicinity of the moving reaction front, 
as well as to a progressive extension of this depletion 
zone (which, recall, is defined as the region where the 
concentrations of the ions are significantly smaller than 
their initial values); of course, as long as permitted by 
the boundaries, on which the concentrations of the ions 
are fixed. This is illustrated in Fig. [T^ 

This depletion effect is more pronounced when the po- 
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FIG. 12; Two snapshots of the concentration profiles of the 
ions in the vicinity of the reaction front at two different times. 
The values of the parameters of the system are: La ~ 1 cm. 
Lb ~ 100 cm, ao = 10 M, and bo = 0.1 M, and applied 
electric field U/L = 1 V/m. 



larity of the applied tension is such that it favors the 
reaction than in the case of the opposite polarity. This 
can be easily understood: in an electric field that favors 
the reaction, not only the disappearance of the ions A~ 
and B~^ is favorized through reaction, but also the back- 
ground ions A~^ and B~ are pushed away from the re- 
action zone. Of course, in a field of reverse polarity, the 
nonreacting ions and B^ arc pushed towards each 
other, and the active ions A^ and B^ are pushed away 
from each other (and thus their consumption through re- 
action is reduced; this leads, in the end, to an extinction 
of the reaction). 

Finally, for a polarity of the applied tension that is fa- 
vorable to the reaction, the expansion of the depletion 
zone is more rapid and pronounced for larger tensions. 
This determines the local increase of the electric field, 
and, in the long time, the numerical instability of the 
integration procedure (as discussed in the previous sub- 
section). 

5. The concentration of the reaction product C 

For the regime of constant jo , as discussed in the main 
text, the concentration of the reaction product c(x) in- 
creases linearly with x with a slope proportional to the 
fieldless value cq, as illustrated by Fig. [T^ 

This initial behavior of the production of C is deter- 
mined by and lasts as long as the initial regime of con- 
stancy of the current density jo, and therefore it is the 
only regime on which we focus in our study. However, 
as discussed above in the corresponding subsection, the 
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current density jo may develop other regimes (if allowed 
by the size of the system). Then, of course, these regimes 
will be reflected by the production of C . In particular, 
the transient increase (if any) in jjoj between the initial 
regime of constancy and the regime jo ~ 1 / Vt is reflected 
by a transient increase in the production of C (a "bump" 
in the profile c{x)). The regime jo ^ 1/v^ itself leads to 
a constant deposition of C in the wake of the front; the 
plateau value of c varies with ao ^md 6o (much in the same 
way as co), but is rather insensitive to the applied ten- 
sion U . Finally, when the reaction front 'hits' the right 
border, there is a great increase in the production of C 
at the border. 
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